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We investigate, if the cosmic ray positron fraction, as reported by the HEAT and AMS collaborations, is 
compatible with the annihilation of neutralinos in the supergravity inspired Constrained Minimal Supersymmetric 
Model (CMSSM), thus complementing previous investigations, which did not consider constraints from unification, 
electroweak symmetry breaking and the present Higgs limits at LEP. We scan over the complete SUSY parameter 
space and find that in the acceptable regions the neutralino annihilation into bb quark pairs is the dominant 
channel and improves the fit to the experimental positron fraction data considerably compared to a fit with 
background only. These fits are comparable to the fit for regions of parameter space, where the annihilation into 
pairs dominates. However, these latter regions are ruled out by the present Higgs limit of 114 GeV from 

LEP. 



1. Introduction 

The cosmic ray positron fraction at momenta 
above 7 GeV, as reported by the HEAT collabo- 
ration, is difficult to describe by the background 
only hypothesis jlj. The data at lower momenta 
agree with previous data from the AMS experi- 
ment J|]. A contribution from the annihilation of 
neutralinos, which are the leading candidates to 
explain the cold dark matter in the universe, can 
improve the fits considerably^,^]. The neutrali- 
nos are the Lightest Supersymmetric Particles 
(LSP) in supersymmetric extensions of the Stan- 
dard Model, which are stable, if R-parity is con- 
served. This new multiplicative quantum number 
for the supersymmetric partners of the Standard 
Model (SM) particles is needed to prevent proton 
decay and simultaneously prevents the LSP a) to 
decay into the lighter SM particles and b) can 
only interact with normal matter by producing 
additional supersymmetric particles. The cross 
sections for the latter are typically of the order of 
the weak cross sections, so the LSP is "neutrino- 
like" , i.e. it would form halos around the galaxies 
and consequently, it is an excellent candidate for 
dark matter. 



In addition to being of interest for cosmo- 
logy, supersymmetry solves also many outstand- 
ing problems in particle physics, among them[|): 
1) it provides a unification of the strong and elec- 
troweak forces, thus being a prototype theory 
for a Grand Unified Theory (GUT) 2) it pre- 
dicts spontaneous electroweak symmetry break- 
ing (EWSB) by radiative corrections through the 
heavy top quark, thus providing a relation be- 
tween the GUT scale, the electroweak scale and 
the top mass, which is perfectly fullfilled 3) it 
includes gravity 4) it cancels the quadratic diver- 
gencies in the Higgs mass in the SM 5) the light- 
est Higgs mass can be calculated to be below 125 
GeV in perfect agreement with electroweak preci- 
sion data, which prefer indeed a light Higgs mass. 
Therefore it is interesting to study the positron 
fraction from neutralino annihilation in the re- 
duced region of SUSY parameter space, where 
these constraints are satisfied and compare it with 
the AMS and HEAT data, as will be done in the 
next section. 
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Figure 1. TTie running of the masses between the 
GUT scale and the electroweak scale. 



2. Neutralino Annihilation in the CMSSM 

In the Constrained Minimal Supersymmetric 
Model (CMSSM) with supergravity mediated 
breaking terms all sparticle masses are related by 
the usual GUT scale boundary conditions of a 
common mass rao for the squarks and sleptons 
and a common mass for the gauginos. The 

parameter space, where all low energy constraints 
are satisfied, is most easily determined by a global 
statistical analysis, in which the GUT scale pa- 
rameters]^] are constrained to the low energy data 
by a x 2 minimization. We use the full NLO renor- 
malization group equations || to calculate the low 
energy values of the gauge and Yukawa couplings 
and the one-loop RGE equations for the sparticle 

lr rhe supergravity inspired parameters are the GUT scale 
Mqut, the common gauge coupling olgut, the common 
mass scales rao, for the spin and spin 1/2 sparti- 

cles, the trilinear couplings Aq and Yukawa couplings Yq 
at the GUT scale of the third generation of fermions, the 
Higgs mixing parameter /i and tan j3 , the ratio of vacuum 
expectation values of the neutral Higgs doublets. 



masses with decoupling of the contribution to the 
running of the coupling constants at threshold. 

Fig. |l] shows a typical running of the com- 
mon masses from the GUT scale to low energies. 
The squarks and gluinos get a higher mass than 
the sleptons due to the gluonic contributions of 
the strong interactions. If m 1 / 2 is not strongly 
above rao, the lightest mass is the supersymmet- 
ric partner of U(l) gauge boson, the bino, which 
mixes with the W 3 boson and spin 1/2 higgsinos 
to neutralinos. The low energy value of the LSP 
is roughly 0.4 times its starting value at the GUT 
scale, i.e. 0.4rai/ 2 . From Fig. |l| it is obvious 
that if rax/2 is much larger than rao, the right 
handed stau becomes the LSP, but this cannot 
be a candidate for neutral dark matter. So from 
cosmology we require the LSP to be neutral, i.e. 
rao sufficiently large compared to rai/ 2 . 

Also the evolution of the mass parameters in 
the Higgs potential are shown using the full 1- 
loop contributions from all particles and spar- 
ticles. The large negative contributions to ra 2 
from the top Yukawa coupling drive it negative, 
thus inducing spontaneous electroweak symme- 
try breaking. The main parameter to get correct 
EWSB is the starting value of ra 2 = {jLq + ra^ at 
the GUT scale. EWSB thus determines the value 
of the Higgs mixing parameter /i 2 for a given value 
of rao- It turns out that as long as rao and ra ly / 2 
are of the same order of magnitude, EWSB re- 
quires the value of \i to be larger or of the same 
order of magnitude than rai/ 2 . 

As a consequence, the Higgsinos are heavy com- 
pared to the LSP, so they hardly mix with the 
LSP, which is then practically a pure bino in 
the parameter space of interest, i.e. for LSP 
masses above 100 GeV, since lighter LSP masses 
do not yield positrons in the interesting momen- 
tum range covered by the HEAT experiment. The 
gaugino fraction is defined as {N^ x + Nf 2 )-> where 
the coefficients determine the neutralino mixing 



As 



shown in Fig. || the gaugino fraction is close to 
one for LSP masses above 100 GeV, i.e. rai/ 2 > 
250 GeV. The gaugino fraction is important, since 
the neutralino properties are quite different for a 
pure gaugino or pure Higgsino. 

Furthermore, the large values of ji imply that 
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Figure 2. The allowed region (green) in the 
CMSSM parameter space for tan/3 =35. Note 
that the combination of low energy constraints re- 
quires m 1 / 2 to be above 350 GeV, which corre- 
sponds to a minimum LSP mass of 150 GeV. For 
lower values of tan (3 the constraints are stronger, 
i.e. the excluded region is larger. 
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Figure 3. The gaugino fractions as function of 
mo and mi/ 2 for tan/3 = 35. For lower values of 
tan (3 the gaugino fraction is even closer to one. 



the pseudoscalar Higgs boson is heavy compared 
with the lightest Higgs boson, in which case the 
latter has the properties of the SM Higgs boson. 
The limit on the SM Higgs boson of 114 GeV from 
LEP restricts the parameter tan (3 to be above 
4.3§. 

The GUT scale and Yukawa couplings are de- 
termined from the requirement of gauge unifica- 
tion and the masses of the third generation parti- 
cles, b — r Yukawa coupling unification can be 
imposed, since these fermions are in the same 
multiplet in any larger gauge group containing 
the SM gauge groups as subgroups. After run- 
ning the Yukawa couplings down to low energy, 
one can indeed obtain the correct b-quark - and 
r-mass for a common Yukawa coupling at the 
GUT scale. But these Yukawas influence the top 
quark as well and the correct top quark is only ob- 
tained for tan (3 =1.6 or tan (3 between 30 and 55, 
where the values above 50 correspond to a posi- 
tive sign of the Higgs mixing parameter fi. The 
positive sign is preferred by the positive deviation 
of the anomalous magnetic moment of the muon 
from the SM, assuming that this 1.6<r deviation 



will be confirmed by future measurements. The 
low tan/3 solution is ruled out, since the Higgs 
limit of 114 GeV at the 95% C.L. from LEP 
requires tan/3 to be above 4.3||, as mentioned 
above, so Yukawa unification suggests the high 
tan/3 solution. 

Nevertheless, we scan over tan/3 between 1 
and 50 to study the different annihilation chan- 
nels, which strongly depend on tan/3 . In ad- 
dition we scan over the remaining parameters 
m o> m i/2- The trilinear couplings shift the Higgs 
mass: Ao > lowers the Higgs mass for given val- 
ues of mo and m^, so the largest excluded region 
is obtained for Ao > 0. Positive values of Ao are 
preferred by the b^X s j data, so if one combines 
all constraints (gauge unification, EWSB, b^X s j 
and a M in a single x 2 function, one can obtain 
a 95% C.L. contour for the excluded region, as 
shown in Fig. || for tan/3 =35. For details we 
refer to previous publications @. 

Neutralino annihilation can occur through 
Z- and Higgs exchanges in the s-channel and 
sfermion exchange in the t-channel. This anni- 
hilation in the halo of the galaxies will produce 
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Figure 4. 77ie neutralino annihilation total cross 
section as function of the center of mass momenta 
of the neutralinos for lepton and W+W~ final 
states. Note the p-wave suppression at low mo- 
menta for light fermions. 
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Figure 5. The neutralino annihilation cross 
section for various diagrams as function of the 
center of mass momenta of the neutralinos for 
tan (3 =35. The curve labeled tot is the total cross 
section including interferences, which is larger for 
bottom than for top final states at large tan (3 due 
to the different signs of the interference terms. 



antimatter at high momenta, thus anomalies in 
the spectra of positrons and antiprotons provide 
an excellent signal for dark matter annihilation. 
The neutralino is a spin 1/2 Majorano particle, 
so it obeys the Pauli principle, which implies an 
asymmetric wave function or spins antiparallel 
for annihilation at rest. This results in a p-wave 
amplitude to the fermion-antifermion final states, 
which is proportional to the mass of the fermion 
in the final state. Therefore heavy final states 
are enhanced at low momenta, as demonstrated 
in Fig. Q for leptons and W-bosons. The same is 
true for quark final states. 

The total cross section is determined by the 
sum of the individual amplitudes; their relative 
signs determine the interferences. The t-channel 
has for t-quark and b-quark final states similar 
amplitudes, but the Z-exchange is for t-quarks 
much stronger than for b-quarks (amplitude oc 
mass). For Higgs exchange the b-quark final state 
becomes large at large tan (3 because of the cou- 
pling proportional to sin/3, while the coupling 
to t-quark final states, proportional to cos/3, is 
strongly suppressed at large tan (3 . Since the Z- 
exchange and Higgs exchange have an opposite 
sign, the t-quark final state is plagued by neg- 
ative interferences, while the b-quark final state 
is enhanced by positive interferences and domi- 
nates at large tan/3 . The various contributions 
are shown in Fig. ||. 

The total annihilation cross section is a strong 
function of the SUSY mass scale, since the t- 
channel contribution is proportional to one over 
the sfermion mass squared, as demonstrated in 
the mo,mi/ 2 plane in Fig. ||. The LSP mass 
is approximately OAmi/2 and for tan/3 =1.6 one 
observes the threshold peaks for 66-, and 
ti final states along the rai/2 axis. For values of 
tan/3 > 5 the bb final state dominates and the 
annihilation cross sections quickly increases with 
tan/3 , as demonstrated by the different vertical 
scales in the figure: at large tan/3 the cross sec- 
tions are several orders of magnitude larger and 
keeps growing with approximately tan 2 /3. The 
boost factors needed for the best fit to describe 
the HEAT data at large momenta are correspond- 
ingly lower, as shown in Fig. 0. These boost fac- 
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Figure 6. The thermally averaged cross section 
times velocity for neutralino annihilation as func- 
tion of mo and for three different values of 
tan/3 (1.6, 5 and 35 from top.) 
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Figure 7. The boost factors from the best fit to 
the combined A MS and HEAT data as function 
of mo and rai/2 for three different values o/tan/3 
(1.6, 5 and 35 from top). 
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Figure 9. The preferred region (light, green) of 
relic density between 0. 1 and 0. 5 for tan (3 =35, 
as calculated with DarkSUSY. The excluded re- 
gion, where the stau would be the LSP, is also 
indicated. For lower value of tan f3 the green re- 
gion rapidly shrinks. Coannihilation of neutrali- 
nos and staus are not considered in DarkSUSY, 
which would extend the green region in a narrow 
strip along the u excl. LSP" region. 
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Figure 8. Fits to the data for tan/3 =1.6 and 
LSP masses of 180 and 130 GeV (2 plots at the 
top) have tt and W + W~ as dominant annihila- 
tion channels, but the fit for large tan /3 with bb 
as dominant channel yields a similar x 2 • 



Figure 10. The x 2 distribution of the fit to pub- 
lished AMS$ and HEATj$ data for tan/5 =35. 
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Figure 11. Expected statistics after one year of 
data taking with A MS- 02 for a neutralino of 240 
GeV and assuming a boost factor given by the best 
fit to the present data. 
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Figure 12. The x 2 distribution of the fit to the 
statistics of the future AMS- 2 data of Fig. [7J 



tors were used in the fit as an arbitrary normal- 
ization, since the dark matter is not expected to 
be homogeneous, but shows some dumpiness due 
to gravitational interactions. Since the annihila- 
tion rate is proportional to the square of the dark 
matter density, the dumpiness can enhance the 
annihilation rate by several orders of magnitude. 

The cross section calculations were done with 
CalcHEP§ and DarkSUSY§, which were found 
to agree within a factor of two for the it final state 
and considerably better for other final states. 

Fig. ^ shows the fit to the data for different 
regions of parameter space with different main 
annihilation channels. The fits were done with 
DarkSusy and follow the same principle as the 
ones from Ref. ||, i.e. the background is left free 
within a normalization error of about 20% and 
the signal can be enhanced by a boost factor as 
discussed above. 

However, we impose the CMSSM constraints, 
which imply tan (3 > 4.3, in which case the bb fi- 
nal state dominates. The positron spectrum from 
b-quark decays can fit the present data as well as 
the harder decays from W-pairs, if the neutrali- 
nos are somewhat heavier, as shown by Fig. ||: 
the fit with b-quark final states yields a y 2 /d.o.f 
of 26/16, which is not much worse than the one 
for W+W~ final states (24/16). These values 
are considerably better than the background only 
fit, which yields 48/16. Here we used the back- 
ground from Moskalenko and Strong as im- 
plemented in DarkSUSY. Since the cross sections 
for bb final states are at least an order of magni- 
tude above the cross sections for ti and at 
large tan f3 , the needed boost factor for the best 
fit is correspondingly lower, especially for lower 
values of mo, as shown before in Fig. [?]. The re- 
gion with large cross section and correspondingly 
lower boost is also the region with a relic den- 
sity parameter between 0.1 and 0.5 of the critical 
density, as shown in Fig. ^. A range between 
0.1 and 0.3 is preferred by the determination of 
the cosmological parameters from the red shift of 
distant supernovae and the acoustic peak in the 
microwave background [11|. 

The fits were repeated for all values of mo and 
rai/2- The resulting x 2 values are plotted in Fig. 
[To[ One observes a fast decrease in \ 2 f° r va l- 
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ues of mi J 2 above 230 GeV, i.e. for LSP masses 
above ~ 100 GeV. Unfortunately the data are 
not precise enough to prefer a certain value of 
the LSP mass. However,the HEAT ballon exper- 
iments correspond to only a few days of data tak- 
ing. With future experiments, like PAMELA(i|] 
on a russian satellite or AMS-02[||| on the ISS 
(International Space Station), one will take data 
for several years, thus being able to decide if the 
excess in the HEAT data above the present best 
background estimate is due to a bad knowledge 
of the background or if it is really a signal for new 
physics. 

If we assume the excess is due to neutralino 
annihilation, then the statistical significance after 
one year of data taking with a large acceptance in- 
strument, like AMS-02|l3|] with an acceptance of 
more than 0.04 m 2 sr, will look like the points in 
Fig. [0]. Here we assumed a boost factor as given 
by the fit to the AMS and HEAT data (Fig. M). 
The resulting y 2 distribution is shown in Fig. [l^. 
Clearly, if the background estimates will be con- 
firmed by accurate measurements of the electron 
spectrum, as will be done by the future expere- 
ments, then the positron spectra can give a clear 
indication of neutralino annihilation with a rather 
precise determination of the neutralino mass. 

We like to thank Drs. L. Bergstrom, J. Edsjo 
and P. Ullio for helpful discussions. 
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